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Abstract —We consider energy-efficient adaptive power al¬ 
location for three incremental multiple-input multiple-output 
(IMIMO) systems employing ARQ, hybrid ARQ (HARQ) with 
Chase combining (CC), and HARQ with incremental redundancy 
(IR), to minimize their rate-outage probability (or equivalently 
packet drop rate) under a constraint on average energy con¬ 
sumption per data packet. We first provide the rate-outage 
probability expressions for the three IMIMO systems, and then 
propose methods to convert them into a tractable form and 
formulate the corresponding non-convex optimization problems 
that can be solved by an interior-point algorithm for finding 
a local optimum. To further reduce the solution complexity, 
using an asymptotically equivalent approximation of the rate- 
outage probability expressions, we approximate the non-convex 
optimization problems as a unified geometric programming 
problem (GPP), for which we derive the closed-form solution. 
Illustrative results indicate that the proposed power allocation 
(PPA) offers significant gains in energy savings as compared to 
the equal-power allocation (EPA), and the simple closed-form 
GPP solution can provide closer performance to the exact method 
at lower values of rate-outage probability, for the three IMIMO 
systems. 

Index Terms —Incremental MIMO, low-complexity MIMO, 
ARQ, HARQ, Chase combining, incremental redundancy, power 
allocation. 

I. Introduction 

Multiple-input multiple-output (MIMO) transmission 
schemes are most suitable for systems with high spectral 
efficiency requirement. Despite having many advantages, one 
of the fundamental limitations of MIMO systems is the cost, 
increased power/energy consumption and the complexity 
associated with their implementation in practical systems 
JT|. Towards addressing these problems associated with the 
conventional MIMO systems, spatial modulation (SM) has 
been proposed in El as a low-complexity MIMO transmission 
scheme that can improve the energy efficiency (EE) with only 
channel state information known at the receiver 0. 

Incremental MIMO (IMIMO) 01, 0 is a variation of SM, 
in which the multiple antennas at the transmitter are used in an 
incremental fashion by utilizing the ARQ feedback to improve 
the reliability. In an IMIMO system, the encoder functionality 
is simplified by letting only one antenna from the transmit 
antenna array to be used to transmit the information at any 
given time. Because of this, only a single RF chain and a 
single power amplifier can be used on the transmitter side and 
the receiver with multiple antennas can decode the message 
optimally with relatively low complexity. After sending the 
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information from a chosen transmit antenna, the transmitter 
waits for the ARQ feedback. If the information is received 
successfully, the receiver sends a positive acknowledgment 
(ACK) and the next packet in the queue is transmitted in the 
next transmission round. If the transmission is not successful, a 
negative ACK (NACK) is sent from the receiver, the same mes¬ 
sage is encoded and sent through a different transmit antenna 
to exploit the spatial diversity. There are three possible ways in 
which the encoding and decoding operations can be performed 
during the transmission of an erroneous packet and they are 
IMIMO using ARQ, CC-HARQ and IR-HARQ, respectively. 
Readers are encouraged to refer to 01, 0 for more details of 
the three IMIMO systems considered in this work and their 
advantages over the conventional MIMO systems. 

Related Work : Previous works on MIMO with ARQ con¬ 
sidered different aspects of the system performance. In 0, 
diversity-multiplexing-delay tradeoff of MIMO ARQ systems 
has been studied. A multi-bit feedback scheme for MIMO IR- 
HARQ was proposed and an outage analysis was presented 
in 0. A progressive ARQ precoder design for MIMO trans¬ 
mission systems to minimize the mean-square error has been 
proposed in 0. The idea of using ARQ feedback for low- 
complexity MIMO system implementation was proposed in 
01, 0 along with an outage analysis of IMIMO systems 
employing three retransmission mechanisms. In 0, among 
other things, the authors showed that for many MIMO-ARQ 
schemes, the efficiency of ARQ protocols is dependent on the 
considered scheme through the accumulated mutual informa¬ 
tion and is independent of the performance metric. 

Recently many works have been focusing on the optimiza¬ 
tion of resources in HARQ systems when the channel state 
information (CSI) is not available at the transmitter. A fixed 
outage probability analysis of HARQ in block-fading channels 
with statistical CSI at the transmitter was presented in 001 . 
Optimal power allocation for improving the average rate per¬ 
formance of HARQ schemes was presented in [ll]] for quasi¬ 
static fading channels with different forms of CSI feedback. 
Power adaptation to minimize the average transmission power 
under a fixed rate-outage probability constraint for both the IR- 
and CC-HARQ schemes was studied in IT2i . A rate allocation 
and adaptation policy based on dynamic programming (DP) 
was proposed in H~3j for truncated IR-HARQ systems. In the 
works of da, mi, the authors proposed power adaptation 
for IR- and CC-HARQ systems in single-input single-output 
(SISO) i.i.d. Rayleigh fading channels to minimize the rate- 
outage probability under an average energy constraint. 

Contributions: We consider the problem of minimizing the 
rate-outage probability of IMIMO systems employing ARQ, 
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Figure 1. System model for the IMIMO. 


CC-HARQ and IR-HARQ under a constraint on average en¬ 
ergy consumption per packet. In particular, we first generalize 
the system model of [j4j to allow for adaptation of transmission 
power in different ARQ rounds and provide the expressions for 
the rate-outage probability of IMIMO employing ARQ, CC- 
HARQ and IR-HARQ. We then formulate the optimization 
problems for each of the three IMIMO schemes using the 
derived outage probability expressions. However, the given 
rate-outage probability expressions are not mathematically 
tractable to be used in an optimization problem formulation. 
Hence, we propose methods to convert these expressions into 
a tractable form and formulate a non-convex optimization 
problem that can be solved by an interior-point algorithm 
for finding a local optimum. To further reduce the solution 
complexity, we propose an asymptotically equivalent approx¬ 
imation of the derived rate-outage probability expressions 
to approximate the non-convex optimization problems as a 
unified geometric programming problem (GPP), for which the 
closed-form solution is derived. 

Even though we consider the same optimization problem as 
in m-tm the present work differs in terms of the system 
model in the sense that here we consider low-complexity 
IMIMO systems which utilize the ARQ feedback to exploit 
the spatial diversity, whereas, in fl4]| - liT6l , point-to-point SISO 
systems with IR and CC-HARQ were considered. 

II. System Model and Rate-Outage Analysis 

We consider a point-to-point IMIMO system having M 
antennas at the transmitter and N antennas at the receiver 
as shown in Fig. [j] We assume a frequency-flat Rayleigh 
block-fading channel. The fading coefficient between the mth 
transmitting antenna and the nth receiving antenna . 1 < 
n < N and 1 < m < M is i.i.d. with distribution CM (0,1). It 
is assumed that h nrn remains unchanged during a fading block 
of a fixed number of transmissions, and change independently 
from one block to another. As in ii, a, we define one 
IMIMO round as up to L possible transmissions for each data 
packet. If the destination is not able to decode a packet after L 
transmission attempts, the packet is dropped. The transmitter 
is assumed to have only statistical knowledge of the fading 
coefficients, whereas, the receiver is assumed to know the 
fading coefficients perfectly. Moreover, in IMIMO systems, a 
new transmit antenna is used for sending an erroneous packet 
in each transmission round (ARQ round). Hence the effective 
channel changes independently within each ARQ round of 


IMIMO and the instantaneous CSI feedback from the receiver 
is not useful. We also assume that L < M so that the quasi¬ 
static Rayleigh block-fading IMIMO system model described 
above can be seen as a single-input multiple-output (SIMO) 
system using ARQ or HARQ, in which the channel fading 
block is equivalent to one ARQ round. 

Each ARQ round consists of T symbols. We assume that 
the modulation symbols have unit average energy and the same 
power scaling factor Pi is applied to all the T symbols during 
the Zth ARQ round. We write the received signal during the 
Zth channel use of Zth ARQ round using the mth antenna for 
transmission as y i [t] = '/Pih rn x rn j [f] + z i [i], where h m = 
[hi ^ m ,..., ZiAr im ] 3 denotes the channel vector from the mth 
antenna to the N antennas at the receiver. The index of the an¬ 
tenna used for transmission in each ARQ round is assumed to 
be known to the receiver. x rn j [f] £ C denotes the modulation 
symbol from the mth transmit antenna in the Zth channel use of 
the Zth ARQ round, y i [Z] = [yij [t],..., [Z]] T denotes the 

channel output and z i [Z] = [zi.i [t ],..., Zn,i [Z]] T represents 
the noise at the receiver and we assume z Ut i [Z] ~ CM (0,1) for 
n = 1,..., N and Z = 1 ,,L. The codebook construction 
and decoding operations for each of the three IMIMO schemes 
has been described in a, G3- Using similar assumptions 
as in M-M about the codewords, we consider rate-outage 
probability defined as the probability that the instantaneous 
rate is smaller than the target rate as a performance metric. 


A. Rate-Outage Analysis of IMIMO Employing ARQ 

For the case of IMIMO employing ARQ, the receiver only 
uses the information from the current ARQ round to decode 
a message. For a target transmission rate of R bps/Hz, the 
probability of outage after Z ARQ rounds is given by: 
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where 7 (s,x) = f* t s ~ 1 e~ t dt is the normalized lower 
incomplete Gamma function, T (N) is the Gamma function 
and Zk = ( 2 R — l) /P k ,l < k < l. The relation in (flbl) 
uses the fact that ||hfc || 2 ~ y 2 distributed random variable 
with 2N degrees of freedom and whose probability density 
function is given by /|| hfc ||2 (fj) = ^ > 0 . i n 

d. -Pmin = min (Pi, • • • , Pi), and we have written the rate- 
outage probability as the sum of the first term and the higher- 
order terms. 


B. Rate-Outage Analysis of IMIMO Employing CC-HARQ 

In case of IMIMO employing CC-HARQ, the receiver 
combines the information received across different transmis- 
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sion rounds using maximal-ratio-combining (MRC). The rate- 
outage probability after l ARQ rounds can be expressed as: 


IMIMO,CC-HARQ A. p 
Pout.Z ~ rr 
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where a k — Pk ||hfc|| ,1 < k < l has Gamma distribution 
with the shape parameter N and the scale parameter P k . 
The term 0 = Y^k=i ak is a sum °f independent and non- 
identically distributed Gamma random variables. Using the 
results from ED, ED, we can express 0 asfl 
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where the notation (xi) = T (xi + t/i) /T (xi) with (xi) 0 = 
1 and the function d 1 ) (.) in ( 13hi is the confluent Lauricella 
hypergeometric function of l variables ED. In ED, we have 
written the rate-outage probability expression as the sum of 
the first term and the higher-order terms. 


C. Rate-Outage Analysis of IMIMO Employing IR-HARQ 

The rate-outage probability after l ARQ rounds for an 
IMIMO system with IR-HARQ can be expressed as J4), 0: 


p“,r- MARQ = r4zios(i+AiM 2 )<4 

= 2*0, (-iiln2)-0,(0) (4) 

where. 


9i (*) 
Qi (*) 
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with u (t) being the unit step function defined as u (f) = 
1, for t > 0 and u(t) = 0, for t < 0, and the symbol 
* represents the convolution operation. The derivation of 


1 For a detailed derivation of the expressions, readers can refer to GH, 03 
and the references therein. 


PcmU 10 ,IR-HAR Q and qi ( t) is similar to the derivations given 
in ElU We provide the derivation of qi (t), 2 < i < l in 
Appendix A. 

We use the Jensen’s inequality in 0 together with the 
results from Section 1II-B1 to write: 

p ,MMO.m-„AEQ >p r | log ( 1+ l^ p t ||ht|| A < |\ 

r ( iJV + i) 

where Z\=Y ( l ) /P u 1 < i < l and Y (l) = l ( 2 R - l). 



III. Optimization Problems and Solutions 

In this section, we first state the general optimization 
problem and describe methods to solve the problem for each 
of three IMIMO systems considered in this work. We define 
the average transmit energy per packet as: 

L 

El A r, IMIMO,ARQ/HARQ type 

-fravg — -t ^ ^PouM-l 
J=1 


We also define the quantity E avg = E avg /T for mathematical 
tractability. Similar to M. ED, we formulate the general 
optimization problem as: 


IMIMO,ARQ/HARQ type 
mm p nl]t r 
(Pi,P2,...,Pl) ’ 

subject to 0 < Pi, for 1 <1<L, (6) 
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E P IMIMO,ARQ/HARQ type , -p 
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A. Solution for IMIMO employing ARQ 

The rate-outage probability expressions for an IMIMO 
system employing ARQ given in (fTbt and (flcl) involve product 
of integrals and an infinite summation, respectively. Hence, for 
mathematical tractability, and to be used in ©, we approxi¬ 
mate ( | I b| i using the standard Gauss-Legendre approximation 
as: 
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where Wi and f are, respectively, the /th weight and the 
/th zero of the Legendre polynomial of order M\ ll22l eq. 
(25.4.30)]. Note that the accuracy of the approximation in 0 
depends on M\. An arbitrarily accurate approximation can be 
obtained by selecting an appropriate value of M\. In practical 
systems using retransmission schemes with a typical value of 
L = 3, outage probability values in the order up to HP' 1 are 


2 The second term in the equation (18) of m should be —e £ (1 — u (—t)) 
instead of —e~ t u (t — 1). 
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of interest Il24ll . and we observed through numerical result^] 
that Mi = 1024 approximates the outage probability values 
with an approximation error smaller than 10 -6 . Using the 
approximation in (0. we can write the optimization problem 
in 0 for ARQ as: 


min 

(Pi,P 2 ,...,Pz.) 


(2T{N)) L 
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subject to 0 < Pi, for 1 < l < L, 
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The optimization problem in 0 is non-convex and hence 
we are not guaranteed to find the globally optimum solution 
to the problem unless an exhaustive search is performed. 
However, nonlinear optimization techniques can be used to 
find a local optimum of 0. We use an interior-point algorithm 
outlined in |25| which uses either a Newton step or a conjugate 
gradient step using a trust region to find a solution. For each 
feasible point (Pi,..., Pl), we need to perform (M\L — M\) 
function evaluations at the zeros of the Legendre polynomial, 
hence the complexity of finding a solution is high. 


B. Solution for IMIMO employing CC-HARQ 

The expressions given in ( l3al ) and ( l3bl ) are not math¬ 
ematically tractable as functions of optimization variables 
(Pi, • • • , Pl). Hence for mathematical tractability, we approx¬ 
imate (l3at using the standard Gauss-Legendre approximation 
as: 
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In this case also, an arbitrarily accurate approximation can be 
obtained by selecting an appropriate value of M 2 . We observed 
through numerical results that M 2 = 512 approximates the 
outage probability values with an approximation error smaller 
than 10~ 6 . Using 0, the optimization problem in (0 for CC- 
HARQ case can equivalently be written as: 
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3 The actual approximation error depends on the 2Mith derivative of 
fk (t) l22l eq. (25.4.30)]. Nonetheless a reference value can be computed 
numerically by generating many random variables and computing the rate- 
outage probability using nl=i Pr (ll h fcl| 2 < z k)- 
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We used the same interior-point algorithm outlined l25l to find 
a solution for (flOt . 


C. Solution for IMIMO employing IR-HARQ 

Using (0, the optimization problem for the IR-HARQ case 
can be written as: 

min 2 R g L (-R In 2) - g L (0) 

subject to 0 < Pu for 1 < l < L, (11) 

L 

P 1 +yP l (2 R gi-i(-R\n2)-gi- 1 (0)) f P given 
1=2 

where the function gi (t ), 1 < l < L is obtained by the con¬ 
volution of the functions q\ (t ), • • • ,qi ( t) defined in Section 
1II-CI We can express this convolution operation in terms of a 
multiple-integral in l — 1 dimensions. Similar to the techniques 
used in Sections fill- Al and 1III-BI we can approximate the finite 
dimensional integrals as finite sums using the Gauss-Legendre 
approximation, or by applying the method described in j26j 
for two dimensions. These finite summations can then be used 
in the optimization problem of (fill and can be solved using 
interior-point methods. 


IV. GPP Approach and Closed-form Solution 

In this section, we provide approximate expressions for the 
rate-outage probability of IMIMO systems employing ARQ, 
CC-HARQ and IR-HARQ to formulate an unified geometric 
programming problem (GPP), for which the closed-form so¬ 
lution is derived. From the rate-outage probability expressions 
for IMIMO using ARQ, CC-HARQ and IR-HARQ given in 
cni. © and (0, respectively, we neglect the higher-order 
terms and write an asymptotically equivalent approximation 
as: 

IMIMO, ARQ/HARQ type Wl 

Pout,; — ^pN ... pwy 

where 

for IMIMO using ARQ 
for IMIMO using CC-HARQ (13) 
for IMIMO using IR-HARQ 

The motivation for the approximation in (IT2l > is: i) as P m ; n = 
min (Pi, • • • , Pi) —> oo , the O (.) terms in (ITdb . (CTcI i and 0 
go to zero faster than the approximated terms in ( IT2l >. ( IT3] >: 
and ii) the maximum possible diversity order achievable in a 
Rayleigh fading channel for IMIMO system with N receiving 
antennas after l ARQ rounds is IN, which is also achieved 
by the approximations in (IT2l) and (fl3l >. Even though the 
asymptotically equivalent approximations of the rate-outage 
probability for the three IMIMO methods have a similar 
structure, they differ in terms of the coefficient Wi as shown 
in ©. The similarity in the structure of approximated outage 
probability expressions in (IT2l) allows us to approximate the 
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optimization problem in HI as a unified GPP as: 


Subject to 


W L 


1 <1<L, 
Wi-1 


0 < Pi, for 
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In the following theorem, we provide the closed-form solution 
for the problem in (1141 ). 


Theorem 1. The closed-form solution for the problem in ( 1741 ) 
is given by: 


P* 
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E^ en N(N + l)^ 1 
(N + 1) L — 1 

Wj -2 /„* NtV+l 

Wi-i (1 + TV) 1 *- lj 
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Proof: Please see Appendix B. ■ 

As can be seen from (fl5l ). the solution of the GPP approach 
differs for different IMIMO methods through the coefficient 
values Wi , 1 < l < L. 


V. Numerical Illustrations and Discussion 

In this section, we present illustrative examples for a perfor¬ 
mance comparison of the proposed power allocation (PPA) and 
equal power allocation (EPA). For the non-convex optimiza¬ 
tion case, we solve the optimization problems ([8]), ( 1 1 Oi l and 
(HU using the interior-point algorithm presented in J25l ; this 
method is labeled as ‘PPA-exact method’ in the plotso For the 
‘PPA-GPP approach’ we solved the approximated optimization 
problem in HU for IMIMO with ARQ, CC-HARQ and IR- 
HARQ, respectively. For the EPA case, we solved for P using 
the additional constraint Pi = • • ■ = Pl = P, in ©, ( I I (it and 

HU- 

Proposed Power Allocation vs Equal Power Alloca¬ 
tion-. Figure [2] shows a performance comparison of PPA 
and EPA for IMIMO employing ARQ, CC-HARQ and IR- 
HARQ under different system parameter values. We plotted 
Pout.z as a function of b given- Following 

observations can be made from Figs. |2(a)]2(c)1 First, for 
higher values of outage probability, the EPA has a similar 
performance as that of the ‘PPA-exact method’, especially 
when the diversity order is high. The gains offered by the 
‘PPA-exact method’ over EPA are more significant for smaller 
values of rate-outage probability (equivalently for higher av¬ 
erage energy limit). Second, for the case of L = N = 2, 
at a rate-outage probability of 10” 4 5 , the gain for the ‘PPA- 
exact method’ over the the EPA solution is 4 dB, 3.1 dB and 
2.3 dB, for IMIMO with ARQ, CC-HARQ and IR-HARQ, 
respectively. However, the gains reduce as the diversity order 
of the system increases (i.e., as the value of N increases^ 

4 Although we cannot guarantee the global optimum, we have verified that 
the solution offered by the interior-point algorithm matches very closely with 
the optimal solution found by an exhaustive grid search over the possible 
values of P;, 1 < l < L. 

5 Note that one can increase the diversity order by increasing the value of L 
as well. However because of space constraints, we could not show the results 
here. 





Figure 2. Performance comparison of the proposed power allocation with the 
equal power allocation for an IMIMO system employing ARQ, CC-HARQ 
and IR-HARQ. The parameters for the simulation are L = 2, M = 2 and 
R = 2 bps/Hz. 
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Figure 3. Performance comparison of PPA-exact method and EPA for IMIMO 
system employing IR-HARQ with different values of spectral efficiency, and 
M = L = N = 2. 

Third, in general, the closed-form ‘PPA-GPP approach’ pro¬ 
vides higher outage probability than the ‘PPA-exact method’ 
and the performance gap between the two proposed schemes 
get closer as the energy limit E glven increases, especially for 
smaller values of L and N. The reason for the “loss” seen by 
the ‘PPA-GPP approach’ relative to the EPA for higher values 
of rate-outage probability is as follows. The approximation 
error of outage probability expressions is non-negligible for 
smaller values of E g i ven . The approximations become tighter 
(asymptotically equivalent) and the performance of the ‘GPP 
approach’ matches that of the exact method as E g i ven value 
increases. To reduce the loss of the ‘PPA-GPP approach’ 
relative to ‘PPA-exact method’, one method is to use tighter 
approximations by considering the higher-order terms of the 
rate-outage probability expressions. However, when higher- 
order terms are also considered, they may include both positive 
and negative terms in the approximations, and this may restrict 
the use of the geometric programming approach to find a 
solution. 

Comparison for Different Values of R: For an IMIMO sys¬ 
tem employing ARQ and CC-HARQ, we have the following 
proposition. 

Proposition 1. In an IMIMO system employing ARQ and 
CC-HARQ, for a given maximum number of transmis¬ 
sions L and target rate-outage probability value of p, if 
(Pi,Ri i P‘ 2 ,R.i j''' ; Pl, Ri ) is the optimal power allocation so¬ 
lution with the average energy E avgt ifor a spectral efficiency 
R,\ > 0 , then the optimal power allocation solutions and the 
average energy for a spectral efficiency of i ?2 f Ri are given 
by 

P1.R2 = Pi,r i y 2Rl __ , l < l < L, ( 16 ) 

E -F ( 2R2 - l \ 

riavg,2 — -C'avg.l l 2 fl i - 1 / 




Figure 4. Comparison of Pi , 1 < l < L values to be used in different ARQ 
rounds for the three IMIMO schemes for a given -Egiven- 

Proof: The proof follows the same arguments as in the 
proof of Proposition 1 in fl31 . ■ 

Hence, for an IMIMO system employing ARQ and CC- 
HARQ, it is sufficient to solve the optimization problem in © 
for a single value of R and scale the resulting power values 
according to (IT6l > to obtain the optimal power values for a 
different value of R. In fact, a similar result as in Proposition 
1 applies to EPA and GPP approaches as well. Hence for a 
given change in the value of R, the performance of all power 
allocation methods shift by the same amount and hence the 
relative performance difference remains the same independent 
of the value of R. However, for an IMIMO system with IR- 
HARQ, optimal power values for different values of R does 
not scale according to (fl6l> . and hence performance difference 
between ‘PPA-exact method’ and EPA is different for different 
values of R, this can be seen from Fig. 0 

Power Values: Figure Q] shows the power values obtained by 
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solving the optimization problems using different approaches. 
As can be seen from Fig. |4(a)| for a given value of E g i ven , 
and for L = 2, the three IMIMO systems have similar 
optimal power values obtained using the exact method, and 
the same power values obtained using the GPP approach. 
The reason for this can be explained by noting that, all the 
three optimization problems have the same average energy 
constraint as p out 1 ' is the same tor all the 

three IMIMO systems (both the exact and approximated 
value). The approximated expression for packet drop prob¬ 
ability pIMIMO,ARQ/HARQ type ^ y (p2p2) for a jj the three 

IMIMO systems, and hence they have the same solution using 
the GPP approach and similar optimal power values using the 
exact method. For L = 3 in Fig. |4(b)| we can clearly see the 
difference in the values of Pi for the three IMIMO schemes 
using the exact method. Furthermore, as the value of E g i ven 
increases, more power is allocated for “later ARQ rounds”, 
which are towards the end of the ARQ process. Since the 
objective of the optimization is to minimize the outage after L 
ARQ rounds, we need to improve the probability of successful 
decoding during these later ARQ rounds. In other words, the 
energy cost associated with an unsuccessful decoding during 
these later ARQ rounds increases. Hence, for large values of 
EgWen, PPA assigns high power values to these later ARQ 
rounds. We can also note from Fig. |4(b)| that for smaller values 
of E g i ven , in case of IMIMO with ARQ, it is optimal to use the 
total transmission power during the first transmission attempt. 

Practical Aspects : For limited real-time computational re¬ 
sources, one can solve the optimization problems offline by 
using nonlinear optimization techniques and store the results 
in a lookup table. For real-time online power allocation, using 
nonlinear optimization techniques to solve © may be too 
costly. In such cases, to achieve low target rate-outage prob¬ 
ability, the simple closed-form ‘PPA-GPP approach’ is more 
computationally efficient and can provide closer performance 
to the ‘PPA-exact method’. If high outage probability (e.g., 
around 10~ 3 or higher) is acceptable, then one can use the 
EPA method. 


VI. Conclusions 

We considered the problem of energy-efficient adaptive 
power allocation in IMIMO systems. In general, these opti¬ 
mization are difficult to solve as the rate-outage probability 
expressions are not mathematically tractable. We developed 
methods to convert the rate-outage probability expressions into 
a tractable form and solved the optimization problems using 
interior-point algorithms. We used asymptotically equivalent 
expressions of rate-outage probability and presented an unified 
geometric programming formulation for which the closed- 
form solution is derived. Possible extensions to the current 
work include, i) considering IMIMO systems with a subset of 
antennas transmitting (as in generalized spatial modulation) 
in each ARQ round instead of a single transmit antenna; 
ii) solving the optimization problems with the objective of 
minimizing the average delay or maximizing the long-term 
average throughput. 


Appendix A 

Proof of (t), 2 < i < l 
Define Qi (s) as in j4): 


Qi (s) = 
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-dr,i = 2,--- ,Z (17) 


where U (a, 6, c) is the Tricomi confluent hypergeometric 
function |23|, and g, (f), i = 2, ■ ■ • ,1 can be obtained as: 
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In ®, we used the relation C 1 ^(r + l) s ^ 

FTi^- 1 ((r + I)*) = (e sln(r+1) ) " 

(t + In (r + 1)) and 5 (.) is the Dirac delta function. 


Appendix B 
Proof of Theorem 1 
We write the Lagrangian function of (IT4l) as: 

W L 


^{Pi,--- , P L , A, fti, ■ • • ,p L ) = 

/ L 
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1=2 M i - 4-1/ / 4=1 

where A,/ti, • • • , /i/, are the Lagrangian coefficients. Since 
the Karush-Khun-Tucker (KKT) conditions are necessary for 
an optimal solution, we have: 

<9Jz? 
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= 0, for 1 <1<L (20a) 
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p*P* =0, for 1 < l < L 
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(20c) 


Since the objective is to minimize Wl/ (Pf ■ ■ ■ P^), 
from (120cl i. it is clear that p* = 0, for 1 < l < L. 
Considering (120ab for l = L and simplifying, we have 
A* = NWl/ (Wl-i {Pl) N+1 ^J- N° w considering 
for 1 < l < L — 1, together with A*, we obtain 
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Substituting for (Pj* +1 ) N+1 in (OH and simplifying, we obtain 
the recursive relation for 2 < l < L as: 


>* _ Wi — 2 /p* \tV+l 

1 ~ (1 + N)Wi-! 1 l ~ l] 


( 22 ) 


Now to solve for P*, we use the fact that the average energy 
constraint should be satisfied with equality at the optimal 
solution, i.e.. 


p;+ 


WiP% 


W 2 P£ 

N ' 


Wl-iPI 


(p*y (PiPi) (p*...p£_!) 


N = E given, (23) 


Now using (l22l > recursively in < 123 b with Wq = 1, we obtain 
the solution for P* as in (ITH . 
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